f(R) Gravities a la Brans-Dicke 
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We extend f(R) theories via the addition of a fundamental scalar field. The approach is reminis- 
cent of the dilaton field of string theory and the Brans-Dicke model. f(R) theories attracted much 
attention recently in view of their potential to explain the acceleration of the universe. Extend- 
ing f(R) models to theories with scalars can be motivated from the low energy effective action of 
string theory. There, a fundamental scalar (the dilaton), has a non-minimal coupling to the Ricci 
scalar. Furthermore beyond tree level actions will contain terms having higher (or lower) powers of 
R compared to the canonical Einstein-Hilbert term. Theories with f(R) will contain an extra scalar 
degree on top of the ad-hoc dilaton and mixing of these two modes around a stable solution is a 
concern. In this work we show that no mixing condition mandates the form Vi((j>)f(R) + V2{<f))R 2 
for the action. 



I. INTRODUCTION 

Recent observations shows that the universe expands 
in an accelerated fashion. This can most easily be ob- 
tained by adding a cosmological constant (or vacuum en- 
ergy). However unfortunately vacuum energy related to 
the present acceleration and the one predicted by theory 
differ by large orders of magnitude. Recently , theories 
of gravity which modify the Einstein-Hilbert action by 
adding terms with different powers of the Ricci scalar at- 
tracted much attention |l| - [28(. The motivation is that 
there is no ab-initio vacuum energy and one is free of the 
mentioned fine tuning problem. 

Another modification of Einstein-Hilbert action is the 
Brans-Dicke model. In that approach there is a new 
fundamental field which couples non-minimally to the 
Ricci scalar and the coupling will yield Newton's constant 
along with the Einstein-Hilbert action if the field stabi- 
lizes at some value. This possible modification may seem 
arbitrary within the Brans-Dicke approach but string 
theory predicts a fundamental scalar called the dilaton 
and its role is very similar. In fact the two models can be 
related by a conformal transformation. So the motivation 
for this approach can be justified within a string theory 
inspired model. Furthermore it is conceivable that string 
theory will predict certain modification in the low energy 
limit of the theory. These modifications will generally be 
coming from beyond tree level calculations and involve 
higher powers of curvature quantities. 

In this work we confine the study to 4 dimensions even 
though our conclusion is shown to trivially generalize to 
d-dimensions. We further consider models which con- 
tains only algebraic functions of the Ricci scalar R. The 
motivation for this is that such models are free of the 
Ostrogradski [2^| (and [23| for a nice review) instability 
infesting theories with higher derivatives. 

In a f(R) theory there will be a scalar excitation. This 



excitation can be viewed to be the df/dR field. Adding 
another scalar on top of this brings in the issue of mixing 
if one expands around a solution with constant fields. As 
we will elaborate in the text, mixing is undesirable for a 
number of reasons so we confine this work to model with 
no mixing around a stable solution for the field equa- 
tions. We show that the no-mixing condition strongly 
constrains the possible models. The general conclusion 
we advocate is that the most general model without mix- 
ing is the one without a pure potential term for the dila- 
ton and R. 

f(R) theories with a an extra fundamental scalar field 
can be shown to be analogous to a bi-scalar tensor grav- 
ity theory via a conformal transformation of the original 
metric (see for example [3(| and [3l[). However, this 
transformation depends on the form of the function / 
and constrains it a little bit. Furthermore in the con- 
formally transformed equivalent tensor bi-scalar model 
the price for equivalency is an accompanying change in 
the original equations of motion for test bodies. This is 
expected since the test particles fundamentally feel the 
original metric. Such an approach based on a conformally 
transformed metric will bring in a theory with an equa- 
tion of state for the scalars that depends on f(R), as for 
instance shown in [32j for pure f(R) theories. There it 
was argued that data capable of selecting the true phys- 
ical frame can be used to discriminate such equivalent 
approaches. Due to these considerations we choose to 
work with the original metric ansatz without resorting 
to the conformally equivalent description. 



II. f(R) GRAVITIES WITH A SCALAR 
COUPLED NON-MINIMALLY 



We choose the following action for this work, 



S = / dx 4 y/—g 



(1) 
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The field equations resulting from this action can be 
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cast as follows, 

V 2 + / = (2a) 

\g^{V<t>? -\v M 4>v v 4> + G MV = o (2b) 

Here subscript <fi (or R) describe the partial derivative of 
/ with respect to <j> (or R) . We have defined the modified 
Einstein tensor as 



Gfiv = JrR^u - ^fg^ + (s M „V 2 - V M V„) f R (3) 

As usual the 0-0 part of the tensor equation can actually 
be traded for the trace equation 

^(V0) 2 + f R R-2f + 3V 2 f R = O. (4) 

A. Constant cj> and R solutions 

In this work we are interested in solutions with a con- 
stant (j> — (j>° and R — R° . This can be achieved if one 
has 



4 = 0, (5a) 
f° R R° - 2f° = . (5b) 

One also has to satisfy the tensor equation 

f o R R o ^-9%lf o = ■ (6) 

Here and the above, the superscripts o in / and its deriva- 
tives mean that they are evaluated at <jf and R° . 

In a cosmological scenario one usually takes the fol- 
lowing for the metric (where we assumed that observed 
dimensions have no intrinsic curvature) 

ds 2 = -di 2 + ]T e 2A *M (da-ija . (7) 

i=l 

In this ansatz the Ricci tensor is diagonal and the ten- 
sor equation written in the orthonormal frame for the 
solution we seek becomes 

f R R°u + \f° = (8a) 

mi - \r = o (8b) 

The meaning of this equation is clear. Unless f° = 
and /jj = the space components of the Ricci scalar 



must be " numerically" the same if they are in general al- 
lowed to be different. It is rather interesting because 
unlike Einstein gravity this seem to imply non-trivial 
isotropization. Thus the solution we are interested in 
becomes R° = g° [/ R°/4 : . From now on we are taking 
Mt)=A(t). 

The solution R° = g°„i?°/4 clearly does not depend 
on the explicit form of the metric ansatz and is quite 
general in f(R) theories. For instance one of the motiva- 
tions to study f(R) models is the existance of de Sitter or 
anti-de Sitter black holes with arbitrary horizon topology 
[H|]- UK- Such solutions must exist even with our inclu- 
sion an extra scalar to the theory. However, since our 
main motivation in this work is the FRW metric ansatz 
(JT]) we do not pursue this avenue further. 

B. Stability of the solution 

If one can find a solution satisfying the conditions one 
has to study whether small perturbations around the so- 
lution are not growing. That is the solution is stable, at 
least in the linear sense. There are two aspects to the 
study of stability. First the metric scale factors Ai(t) 
must not grow under arbitrary perturbations around the 
solution A°(t). Second the two scalars that are in the 
theory must not have negative mass squared. One of 
these scalars is of course (j>, the other is the scalar that is 
activated in f{R) theories, namely f R . 

Under small perturbations we have 

and 

3f° RR V 2 5R + - (f° R - f° RR R°)5R - f° R4 H . 

This can be cast into the following form 

VHX = SSX ^ 

J RR 

with 5X T = (5R, 8(f)) and the elements of stabilization 
matrix S is given by 

S RR = f R -f° RR R° + 3f R * 2 (10a) 

= -m*f°RR ( 10b ) 

S R * = f° m [-if + 3/^] (10c) 
S* R = -3f; R f° RR (10d) 

For the stability of the <f> and f R excitations the matrix 
S must have positive singular values regarding the fact 
that it is in general not symmetric. 

Apart from requiring positive eigenvalues for S there 
is also the issue of mixing. The field <p is an analogue of 
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the dilaton field predicted by string theory [35|. Dilaton 
stability is an important aspect of string theory moti- 
vated models. For example a rolling dilaton would imply 
a time-dependent Newton's constant which is subject to 
stringent bounds. These bounds can be met within a 
specific model, however for purely acsthctical reasons we 
would like to seek a possibility where it never occurs. On 
the other hand if there is a mixing between the dilaton 
and fn around a constant <j> solution this would imply 
that taking the dilaton to be a constant at the action 
level and assuming that the resulting theory is a stable 
theory would be contradicting arguments. Our main em- 
phasis in this discussion is that the field <fr is supposed 
to be a fundamental scalar field as opposed to the fact 
that fpt excitations are genuinely gravitational in origin. 
We would like to stress this point further by invoking 
the standard model of particle physics. Gravity is as- 
sumed to be a singlet under the electro-weak symmetry 
group and this fact is also consistent with the anomaly 
cancellations. However if one would like to contemplate 
the possibility that (f> is not a singlet under the electro- 
weak symmetries, a mixing between (f> and fa can be very 
problematic. 

So it is somewhat desirable to prevent mixing. As it is, 
there seems to be no natural mechanism [36| to prevent 
4> and fn mixing unless /j^ = f£ R = 0. This can be 
trivially achieved if there are no non-minimal couplings 
between <j> and R but this will immediately take us out- 
side the emphasis of this work. Another instance of no 
mixing can be provided if the extremum values <f> and 
R a are not related by the conditions but this can trivial- 
ize the model. We will see in the next sections that there 
is a way to have no mixing without fine tuning. 

Another important point is that the excitations 5(f> and 
SR must not have vanishing masses as this is strongly 
disfavored from a phenomenological point of view. 

There are various terms that play a privileged role in 
finding solutions, providing stability and preventing mix- 
ing. Terms that depend only on R can have no bearing on 
the extremum condition in the (j) equations. Pure <f> terms 
on the other hand would be important for both sectors. 
Pure <j> and pure R terms cannot directly influence the 
easiest no-mixing condition f^, = 0. Furthermore the 
only term that would not affect the extremum condition 
of the trace equation and the cj> equation is the [globally) 
conformally invariant R 2 term. This term will only influ- 
ence the stability condition. As such it is an important 
term to possibly save a model without stability as the 
role of the R 2 term to save R + 1/R gravity models. 



III. MODELS INVOLVING R p WITH p < 1 

In this section we consider several models of interest 
involving powers of R such that they are more important 
near R = 0. These models attracted much attention re- 
cently in view of their potential of explaining the present 
acceleration of the universe. 



A. Models That Won't Work 

As a first try we study the following 

/(</», R) = a4> 2nl+2 R - (34> 2n2+2 R- m + 7 i? 2 (11) 

Since the R term is not present on its own this is a 
flavor of the Brans-Dicke approach. We have also tried 
to keep the 4> — > — cj> symmetry of the kinetic term of the 
scalar field. It is special in the sense that it does not 
have a pure <fi potential and hence no analogue of a vac- 
uum energy (cosmological constant) . This is aesthetically 
favorable since the main purpose of f(R) theories is to 
explain the acceleration without a cosmological constant. 
As we have mentioned R 2 term has no role in finding an 
extremum and can possibly save the model if without it 
there is instability. The extremum conditions are easily 
met if 



iT'+V" 1-2 " 2 = (2 + to)-, (12a) 

a 

(n 2 + l) = (m + 2)(m + l) . (12b) 

We immediately see one important feature. The ex- 
tremum condition for the trace equation is trying to con- 
formalize the extremum and bringing a constraint on n<i 
and ri\. By conformalization we only want to make a 
reference to the fact that R 2 term cannot influence the 
extremum condition of the trace equation and that it is 
conformally invariant. The extremum conditions are do- 
ing their best to resemble the rest of the action to that 
form. Of course such an effect cannot be present in a 
theory where there is no <fi field. Since n>2 and ri\ are 
related this has direct bearing on the dimensionalities of 
the coupling constants. In fact one can show that for 
all to one has det(S) = [37| meaning there will be a 
massless excitation which is not phenomenologically fa- 
vored. This problem has occurred because there was no 
independent mass scale in the theory. Since there are two 
conditions for extremum one of them served to relate the 
dimensionalities of a and (3. 

We therefor have to add terms which introduces two 
ad-hoc mass scales which will not be constrained by the 
extremum conditions. To exemplify the need for this let 
us try the following where there is only one possible mass 
scale, 

/(</>, R) = ol<\> 2 R + A</> 4 - /3<j> 2n+2 R- ,n + 1 R 2 . (13) 
The extremum conditions give 

-aR- 2\4> 2 + (to + 2)/3(f> 2n R- m = , (14a) 
aR + 2\<p 2 - (n + l){3<j) 2n R- m = 0. (14b) 

which implies 



4 



m + 2 = n + 1 
and again <j> and i? Q cannot be fixed separately. 

B. A Possible Model 

It is now apparent that we have to introduce two mass 
scales that are not related via the extremum conditions. 
If one adds a mass term for the scalar field and the usual 
Einstein-Hilbcrt term this will be trivially done. But we 
would be introducing an analogue of a vacuum energy 
and also go outside the Brans-Dickc approach. One way 
out of this dilemma that might generally work is to never 
allow pure potentials for the cj> field, that is every term in 
the action will have to be an interaction between <j> and 
R. For example the following 



1 



f(<f>,R) = a<f> 2 R- 
could work because one fixes 



>+2^ +a; ^p+2 jR 2 +7jR 2 _ (15) 

R 



R 

/2n+4p 



(n- 2)a 
3w(p+l) ro 
(n - 2) 2 a 3 
27(3to 2 (p+l) 



-2p 



2 ' 



(16a) 
(16b) 



as usual 7 term can be used to save stability. We would 
like to focus on this model a little more. Taking a > 
is important because this term is related to the effective 
Newton's constant. On the other hand, since we would 
like to have R > to account for the acceleration of the 
universe and that 2n + 4p is even we get the following 
bounds on the parameters 



> 0, 



(n-2) 



> 



(17a) 
(17b) 



w(p+l) 

Using the extremum values one can check the mixing 
status of the model. It is rather interesting that 



No Mixing 



n = 



That is the term with the 1/R will have to be a di- 
mension zero operator, the most relevant term possible 
in terms of dimensional analysis. Dimension zero op- 
erators can only occur when there are more than one 
field around and possible inverse powers are in the game. 
Since this is not a fine tuning on the coupling constants 
it is legitimate to assume it. This will imply 



(3 > 
w(p+l) < 



(18a) 
(18b) 



The absence of mixing makes the stability analysis 
much simpler since S is diagonal from which we can read 
the masses as follows, 



JRR 



R° 



m 4> 2 = ~ofu 



(19a) 
(19b) 



The general condition for positivity of the mass of the <j> 
field yields p > 0. The other condition requests that 7 > 
7o where 7 Q depends on the parameters of the system. 
Since p > the positivity of R a imposes to < 0. This 
is only fair because it is the term with the largest power 
of (f) in the model and if we consider the potential for a 
given value of R it must be bounded. 

To conclude the study of this model we cast it into a 
standard form where every parameter is positive and the 
dimensionalities are apparent 



/(</>, R) = a<fR 



R 



P2 



2m+2 



R 2 + cR 2 



Here all the couplings a, b, w and c are dimensionless. 
The dimensions are described by [i\ and \xi. Let us fur- 
ther investigate the behaviour of the model for the par- 
ticular instance of m = 1. In this case we have 



Ro — 



M2-/s , 



m 



fR 



s = 




(20a) 
(20b) 
(20c) 

(20d) 
(20e) 



If one fixed the Hubble constant H from experiments 
one fixes R and hence fii, this number is rather small. 
It follows immediately that one of the excitations will 
have to have a mass of the same scale as the Hubble 
constant. Such a particle if it exists must couple very 
weakly to matter. In fact this somewhat mandates to 
assume that it is a singlet under electroweak and strong 
interactions and therefor couples to the ordinary fields 
only via gravity. This possible handicap actually exists 
in all f(R) theories with or without a Brans-Dicke type 
approach. However in a pure f(R) approach one is not 
forced to assume the excited field fa is a singlet under 
standart model interactions: it is already a singlet from 
the start. So the general conclusion is that the Brans- 
Dicke type field must be a singlet under standard model 
unless rather contrived and ad-hoc hierarchies are im- 
posed. Of course in this discussion we have assumed 
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that the dimensionless coupling a is of the order unity 
we would like to emphasize that if a is allowed to vary 
from unity, as apparent from (|20p . R Q and can be 
very different from each other. This of course will bring 
in the issue of naturalness similar to the u> parameter of 
the standard Brans-Dicke model. 



C. Most General Model Without Mixing 

In an attempt to generalize the model presented above 
one can try the following 

f((t>,R) = acj) 2 R- f3(t) 2n+2 R q -uj<j) 2p+2 R 2 + 1 R 2 . (21) 

with q < 1. A rather interesting feature is that if one 
searches for a no-mixing condition we get the same result, 

No Mixing n = for any q 

The rest of the analysis is similar. The requirement 
that R > will mean li > 0. The stability conditions 
will fix p > and 7 > 7o where 7 is model dependent. 

So if one has a working pure f(R) mode, the quest to 
embed it in a Brans-Dicke approach must start with the 
following 

R) = cf> 2 f(R) - lo^+ 2 R 2 + 7 i? 2 . (22) 

The reason why the method exposed works is of course 
due to the fact that we have used the conformally invari- 
ant term R 2 which will not enter the extremum equations 
on R. So a working f(R) gets us half the way without 
any problem. In the absence of V the condition of no- 
mixing and the extremum equation on cj> are identical so 
we get two birds with one stone. 

In fact the most general model without mixing is the 
following 

f{4>,R) = V x {4>)f{R) + V 2 {4>)R 2 (23) 



IV. CONCLUSION 

In this work we have studied a way to get f(R) gravi- 
ties via a scalar field which could be motivated by the 
dilaton of string theory. We have made emphasis on 
avoiding the mixing between the fundamental scalar and 
the scalar fn that is excited in an f(R) theory. One could 
of course relax this requirement and still can find a viable 
model. However we have shown that the no mixing re- 
quirement yielding the most general model actually kills 
two birds with one stone. That is the extremum condi- 
tion on the fundamental scalar of the model and the no 
mixing condition become identical and hence there is no 
fine tuning of parameters. 

Even though we have confined our study to 4- 
dimensional models generalization of the no-mixing 
model can be easily done. The following will be free of 
mixing in d-dimensions. 



/(</>, R) = Vx(<t>)f(R) + V2{4>)R dj2 (24) 



We have shied away from adding a pure cj> potential on 
grounds that it would complicate the models beyond the 
scope of this work. Furthermore adding such a term also 
opens the question about a pure R term (other than R 2 ) 
and this will take us outside the Brans-Dicke approach. 
Nevertheless it is an important extension of the models 
presented here. 

In this work we have worked with a usual kinetic term 
for the <f> field. It will be interesting to use the canonical 
form of the dilaton action coming from string theory and 
extend the study there. The motivation for this will be to 
provide a stabilization mechanism for the dilaton. So far 
in the literature the majority of approaches are focusing 
on dealing with the stability of the dilaton via energy- 
momentum tensors of compact objects around extra di- 
mensions. It could be interesting to study the impact of 
a f(R) approach in string inspired dilaton theory as it 
could yield non-canonical ways to stabilize the dilaton. 
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Even though we did not take the canonical form of the 
low energy string action our field <p can be considered 
to be directly related to the dilaton of string theory by 
a conformal transformation which is well behaved gen- 
erally in contrast to the possibility of not being able to 
perform the conformal transformation for every / in a 
higher order gravity theory 

That is without contrived cancellations between the 
quantities f°, etc. 

This is also clear from the fact that the extremum con- 
ditions do not fix R° and (f>° separately 



